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Abstract 

\ We simulate the future evolution of the observed inhomogeneities in the local uni- 

J> ■ verse assuming that the global expansion rate is dominated by a cosmological con- 

stant. We find that within two Hubble times (~ 30 billion years) from the present 
\ epoch, large-scale structures will freeze in comoving coordinates and the mass dis- 

tribution of bound objects will stop evolving. The Local Group will get somewhat 
closer to the Virgo cluster in comoving coordinates, but will be pulled away from 
the Virgo in physical coordinates due to the accelerated expansion of the Universe. 
In the distant future there will only be one massive galaxy within our event hori- 
zon, namely the merger product of the Andromeda and the Milky Way galaxies. 
All galaxies that are not gravitationally bound to the Local Group will recede away 
from us and eventually exit from our event horizon. More generally, we identify 
the critical interior overdensity above which a shell of matter around an object will 
remain bound to it at late times. 
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1 Introduction 



Recent data on the temperature anisotropies of the cosmic microwave back- 
ground (e.g. Hanany et al., 2000; de Bernardis et al., 2000), the luminosity 
distance to Type la supernovae (Perlmutter et al., 1998; Riess et al., 1998; 
Garnavich et al, 1998), and the large-scale distribution of galaxies (Peacock 
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et al., 2001; Verde et al., 2002) favors a flat universe with present density 
parameters of fiju.o — 0.30 ± 0.15 in matter and fl\,o — 1 — &>m,o in a cosmo- 
logical constant. The energy density of the vacuum (the so-called cosmo logical 
constant or "dark energy") appears to be currently dominating the expansion 
rate of the Universe. 

Given these specific values for the cosmological density parameters, it has now 
become possible to predict quantitatively the future evolution of the visible 
Universe. The existence of a cosmological constant has profound consequences 
in this context (see, e.g. Starobinsky, 2000; Gudmundsson & Bjornsson, 2002; 
Loeb, 2002); in particular: (i) when the Universe will age by a factor of a few, 
the event horizon will stall at a fixed proper distance of 3.6/ig 1 Gpc around 
us (where ho is the present Hubble constant in units of 100 km s _1 Mpc -1 ); 
(ii) all sources with present redshifts larger than 1.8 have already crossed our 
event horizon and are therefore out of causal contact today; and (iii) even if 
we continue to monitor the sources with current redshifts of z = 5 — 10 into 
the infinite future, we will only be able to see these sources acquire intrinsic 
ages of 4 — 6 billion years in their rest frame. 

In this paper we use N-body simulations to calculate the future gravitational 
growth of the observed large-scale structure in the local universe. For simplic- 
ity, we assume that the vacuum energy density is constant in time (as for a 
classical cosmological constant) and adopt the values Qm,o — 0.3, Qa,q = 0.7 
and h = 0.7. We simulate the evolution of density inhomogeneities within a 
sphere of radius ~ lOO/i^ 1 Mpc around the Milky Way galaxy. Our goal is to 
determine the time when the accelerated expansion of the Universe will freeze 
the large-scale structures in distant future, as well as to find the minimum 
overdensity of matter interior to a shell surrounding an object today that will 
allow the shell to remain bound to the object at late times despite the repul- 
sive gravitational force of the vacuum. Another issue of particular interest is 
the question whether the Local Group of galaxies is bound to the nearby Virgo 
cluster. It has been well-established that the Local Group has a peculiar infall 
velocity towards the direction of the Virgo cluster, the so-called 'Virgocentric 
Infall' (e.g. Aaronson, et al., 1982; Davis & Peebles, 1983; Lynden-Bell, et al., 
1988), but it has not been investigated yet by a direct numerical simulation 
whether the Local Group will eventually fall into the Virgo cluster in the 
future given the observed galaxy distribution in the local universe. The possi- 
bility remains that the Virgocentric Infall will freeze in comoving coordinates 
due to the exponential expansion of the universe, and the Local Group will 
not ultimately merge with the Virgo cluster. We investigate this issue for the 
first time by running an N-body simulation, starting from initial conditions at 
z=0 that match the observed galaxy distribution in the nearby universe and 
ending it in a distant future. 

The outline of this paper is as follows. In § 2 we describe the characteristics of 
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our N-body simulations for the evolution of the global structure of the local 
universe. The results from the simulations are discussed in § 3. We then derive 
in § 4 the present-day overdensity threshold for a spherical region to collapse 
in the future, and apply this condition to the Local Group in § 5. Finally we 
summarize our main conclusions in § 6. 



2 Simulations 

The simulations were carried out with the parallel tree N-body/SPH code 
GADGET 3 (Springel, Yoshida, & White, 2001). The SPH component of the 
code was turned off as we are only interested in the dark matter particles which 
dominate the fluctuating mass density. For the initial conditions, we used the 
simulation output of Mathis et al. (2002) 4 at z — 0. This simulation was de- 
signed so that at z = the simulated mass density field would best match the 
observed galaxy overdensity distribution in the IRAS 1.2 Jy survey of Fisher 
et al. (1994, 1995). Mathis et al. (2002) constructed mock catalogues of the 
Mark III Catalogue of Peculiar Velocities (Willick et al., 1995, 1996, 1997) us- 
ing their simulation result, and have shown that the radial peculiar velocities 
in mock and real catalogues agree very well. The original data set of Mathis et 
al. (2002) contains 50.7 million dark matter particles of mass 3.6 x 10% 1 M Q 
within the comoving radius of 80/1(7 1 Mpc around the supergalactic center (in- 
ner high-resolution region), and 20.5 million particles within the comoving 
radius of ~ 200/iq 1 Mpc (outer low-resolution region). The simulation of this 
spherical region was run with a vacuum boundary condition. 

Because we are only interested in the evolution of the structure with mass- 
scale larger than that of the Local Group, we reduce the original data set by 
a factor of 100 through random sampling of particles. This reduction of data 
allows us to run the simulation with a relatively low computational cost while 
preserving the mass and spatial resolution we need to make predictions for 
the questions we are interested in, given the observed galaxy distribution in 
the nearby universe. The resulting mass of each dark matter particle in the 
high resolution region is 3.6 x l0 11 h () 1 M Cl) (comparable to the mass of an L+ 
galaxy), and the number of high resolution particles is half a million. We have 
set the gravitational softening length to the physical scale of lOO/ig 1 kpc, since 
structures below this scale are not of interest. This mass and spatial resolution 
is sufficient for following the overall structure of the local universe, but not 
sufficient for resolving the inner details of galaxy groups such as the Local 
Group, of which the Milky Way galaxy is the dominant member. 
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Using the reduced data set, we have simulated the evolution of the local uni- 
verse from the present time (corresponding to a scale factor value of a = 1 at 
z — 0) to the 6 Hubble times into the future (a = 166; ~ 84 billion years from 
today). Because the present Hubble time t H = l/H = 1A Gyr (for h = 0.7) 
is very close to the present age of the Universe t = 13.5 Gyr in the adopted 
cosmology, the epoch of a = 166 roughly corresponds to t — to + 6t# ~ 7t# 
from the Big Bang (t — 0). Note that we denote all present-day values with a 
subscript zero. 

In addition to the above constrained realization simulation of the local uni- 
verse, we have also performed a simulation with a periodic box of L = 100/ig 1 Mpc 
containing N = 64 3 dark matter particles, in order to make sure that the sta- 
tistical results presented in this paper are not affected by cosmic variance or 
by the choice of the vacuum boundary condition for the constrained realiza- 
tion run. Statistically, we find that the results from the two runs agree very 
well with one another, and so we only show the results from the constrained 
realization run in § 3. 

The simulations were performed on the local Beowulf PC cluster located at 
Harvard-Smithsonian Center for Astrophysics. It is interesting to note that 
the numerical code speeds up its calculation at late cosmic times. The acceler- 
ated expansion of the Universe freezes the large-scale structures in comoving 
coordinates, and so at late times the variable time step of the code increases 
and the tree update for the force computation becomes less frequent. 



3 Future Evolution of Large-Scale Structure in the Local Universe 

Figure 1 shows the distribution of particles in a slab covering a width (centered 
at zero) along the supergalactic Z axis of —15 < SGZ < 15/iq 1 Mpc, and 
projected onto the supergalactic XY plane in comoving coordinates. From top 
left to bottom right the four panels show snapshots at times t = t , t + t#, 
t + 2tn, and t + 6tn, corresponding to a = 1.0, 2.5, 5.8, and 166, respectively. 
The solid circle in each panel denotes the physical radius of 100/iq 1 Mpc at 
each epoch. At t = to + Qtn, the Universe has expanded so much that this 
circle is no longer visible. Three major density peaks can be easily identified 
in all panels: the Coma cluster at the top, the Great Attractor (Centaurus) 
to the left of the center, and the Perseus-Pisces to the right of the center. 
The Virgo cluster is also apparent just above the center. Close examination of 
the structures reveals that most of the evolution takes place between to and 
to + 2t#, and little evolution is seen after t = to + 2£#. This follows from the 
fact that after a few Hubble times, the exponential growth of the scale factor, 
a(t) oc exp(^l\^t/t H ), damps any peculiar velocity which is gravitationally 
induced by mild density inhomogeneities. The bottom right panel shows our 



4 



event horizon by the thick dashed circle around the supergalactic center, which 
is located at a physical radius of 3.6/iq 1 Gpc. At t = to + 6t#, all clusters 
except Virgo are already outside our event horizon. The Virgo cluster will 
also exit from the event horizon at t « t + Q.3t H (a = 210; ~ 88 billion 
years from today) and its image will fade quickly while remaining frozen on 
the sky at the time of its exit. Although our calculation does not take general 
relativistic effects into account, we expect our near- horizon results to be valid, 
since structures on this scale should at any event be frozen. 

The lack of late evolution is even more apparent in the mass distribution of 
bound objects. Figure 2 shows the cumulative number of dark matter halos 
with mass above M per comoving volume, N(>M). The dark matter halos 
are identified through the HOP grouping algorithm (Eisenstein & Hut, 1998) 
with the outer overdensity threshold of 5 out = 80 and the peak overdensity 
threshold of <5 pcak = 240 in comoving coordinates. From left to right, the lines 
refer to times t = t (solid), t + tjj (short- dashed), t Q + 2£# (dotted), t + 
(long-dashed), to + 4£# (dot-short dashed), to + 5tn (dot-long dashed), and 
to + Qtn (short dash - long dash), which correspond to scale factor values of 
a = 1.0, 2.5, 5.8, 13.5, 31.2, 72.0, and 166, respectively. The last four lines 
are almost identical within sampling errors, showing that the mass function 
freezes after two Hubble times from today and that mergers are very rare 
subsequently. 



4 Gravitationally Bound Objects 

Next we consider the interior overdensity threshold above which a shell of 
matter around a spherically symmetric object will remain bound to it at late 
times. In the absence of a cosmological constant and with no shell crossing, 
the threshold simply amounts to the kinetic+potential energy of the shell 
being negative (Gunn & Gott, 1972). The existence of a cosmological constant 
makes the interior gravitating mass of the vacuum grow larger as the shell 
radius increases; consequently a shell with a negative energy today may still 
be pulled away from the object by the cosmic acceleration at late times. We 
would like to find the minimum mean interior overdensity today, 5 C , above 
which the shell will remain bound to the object at arbitrarily late times. Here 

5 = (p — p)/p, where p = f2 M o(3-£fo/87rG) is the mean density of matter in 
the Universe. 

The critical mean overdensity interior to the shell radius, 5 C , was derived by 
Lokas & Hoffman (2002) in their equation (27). Taking the initial time as 
z = 0, this equation yields 5 C = 17.6 for the set of cosmological parameters 
adopted in this paper. We note that Lokas & Hoffman (2002) have ignored 
the possibility that the mass shell may have a nonzero initial peculiar velocity. 
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Fig. 1. Future evolution of the local universe in comoving coordinates. Shown are 
the particles in a slab of thickness —15 < SGZ < 15/ig 1 Mpc projected onto the 
supergalactic XY plane at the times t = to, to + tu, to + 2i#, and to + Qtn (cor- 
responding to a = 1.0, 2.5, 5.8, and 166) from top left to bottom right. The thick 
solid circle in each panel indicates the physical radius of 100/iq 1 Mpc around the 
supergalactic center. In the bottom right panel, the Universe has expanded so much 
that this circle is no longer visible. Instead, we show the event horizon at a phys- 
ical radius of 3.6/iq -1 Gpc as the thick dashed circle. [The number of particles in 
this figure is further reduced by a factor of 20 from the original version for the as- 
tro-ph submission. See http://cfa-www.harvard.edu/~knagamine/LocalGroup for 
high resolution figures.] 
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Fig. 2. Evolution of the cumulative number of dark matter halos with mass above 
M per comoving volume, N(>M). From left to right, t = to (solid), to + tfj 
(short- dashed), to + 2tn (dotted), to + 3tn (long- dashed), io + 4i# (dot-short dashed), 
to + 5tn (dot-long dashed), and io + 6i/f (short dash - long dash), which correspond 
to a = 1.0, 2.5, 5.8, 13.5, 31.2, 72.0, and 166, respectively. The last four lines are 
almost indistinguishable, implying that the mass function freezes after ~ 2tf{ from 
the present time and the merger or accretion rates are very low subsequently. 

For the growing mode of a density perturbation, the inward peculiar velocity 
will induce collapse at a lower initial overdensity than that found by Lokas & 
Hoffman (2002). 

We have attempted to test the validity of the critical overdensity threshold 
S c = 17.6 with our N-body simulation. Figure 3 shows the physical radial 
velocity of particles in and around the two most massive objects in the sim- 
ulation: Perseus and Centaurus clusters. For the HOP grouping parameter of 
<W = 80 (see § 3), the total grouped masses of Perseus and Centaurus are 
1.6 x 10 15 ho l M Q and 9.0 x 1O 14 /iq 1 M at t = t , and 5.5 x IO^/^Mq and 
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3.0 x 10 15 h 1 M Q at t = t + 6t H , respectively (see § 4 of Mathis et al. (2002) for 
a detailed comparison between the simulated and observed cluster properties). 
On the left column of Figure 3 we show the mean physical radial velocity of 
concentric shells around the density maximum as a function of their physical 
radius at t = t + 6t H . The short-dashed line indicates the Hubble law at 
t — to + 6tff, and the horizontal dotted line indicates zero velocity. Inside the 
radius of 2/iq 1 Mpc from the cluster center the velocity fluctuates around zero, 
implying that the particles have virialized. At larger radii, the velocity field 
approaches the Hubble flow (denoted by the dashed line). Because the phys- 
ical radius of 10/ig 1 Mpc at t = t + Qt H corresponds to comoving 60/ig 1 kpc, 
dark matter particles at that radius are still infalling towards the density peak 
and so the physical radial velocity does not perfectly match the Hubble flow. 

On the right column of Figure 3 we have selected out all particles that are 
within a physical radius of 2h^ 1 Mpc around the density maximum at t = 
to + 6t#; these particles can be regarded as bound and virialized based on 
the left column plots. We then trace back the overdensity that these particles 
had at t — t , and plot it against the radial particle velocity at t — t + 6t#. 
Zero velocity is indicated by the vertical long-dashed line, and the critical 
overdensity 8 C = 17.6 is indicated by the solid horizontal line. For a small 
number of particles, the value of 5 C = 17.6 is an overestimate, probably due 
to them having inward initial peculiar velocities. On average, the analytic 
estimate for the critical overdensity based on the spherical tophat collapse 
model appears to provide a good approximation to the actual threshold. 



5 Future Evolution of the Local Group 

5. 1 Will the Local Group Fall into Virgo ? 

Although our simulation does not resolve the inner structure of the Local 
Group, we can get an idea about the future trajectory that its center-of- 
mass will follow by tracing the particles which reside near the supergalactic 
center at the present time. In Figure 4, we show the projection of a region 
with —15 < SGZ < 15/ig 1 Mpc onto the supergalactic XY plane in a similar 
fashion to Figure 1. The density peak near the top of the panel is the Virgo 
cluster, which is separated by a comoving distance of 14/iq 1 Mpc from the 
supergalactic center. The observational estimates of the distance to the Virgo 
cluster from us by the methods of surface-brightness fluctuations, planetary 
nebula luminosity function, Tully-Fisher relation, and Cepheid variable stars, 
all agree with 16 ± 2 Mpc (e.g. Jacoby et al., 1992; Graham et al., 1999), in 
reasonable agreement with the simulated distance within 2 — a. 
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Fig. 3. Left column: mean physical radial velocity v r of concentric shells around the 
maximum density peak is shown as a function of physical radius at t = to+6tji- The 
short-dashed line indicates the Hubble law at t = to + 6tn- The velocity fluctuations 
around zero inside the radius of 2h^ 1 Mpc from the cluster center indicate that the 
associated particles have virialized. Right column: the present-day overdensity of 
particles that are within the physical radius of 2/iq 1 Mpc around the maximum 
density peak at t = to + 6t# versus the physical radial velocity of these particles 
at t = to + 6iff. Zero velocity is indicated by the vertical long-dashed line, and the 
critical overdensity 5 C = 17.6 is indicated by the solid horizontal line. 

We identified all particles enclosed within a sphere of radius 2/ig 1 Mpc around 
the supergalactic center at the present time as shown in the left panel of 
Figure 4). We have found 11 such particles (shown as the open circles), cor- 
responding to a total mass of 3.9 x 10 12 ho l M & . For comparison, the Local 
Group has a zero-velocity surface (separating it from the Hubble expan- 
sion) of 1.18 ± 0.15 Mpc and a total mass of M LG = (2.3 ± 0.6) x 1O 12 M 
(Courteau & van den Bergh, 1999; van den Bergh, 1999) based on the ra- 
dial velocity dispersion of its members. Initially, the sampled 11 particles 
have a mean peculiar velocity of 670 km s" 1 at t = to towards the direc- 
tion of (—250, 502, —368) km s _1 with respect to the cosmic microwave back- 
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ground (CMB) frame. For comparison, the observed motion of the Local 
Group relative to the CMB frame is 627 ± 22 km s _1 in the direction of 
(—406, 352, —324) km s _1 (Kogut et al., 1993). The agreement between the 
simulated and the observed velocity amplitude at t — t is impressive given 
the crudeness of the initial conditions of the simulation; however, the direction 
of the velocity vector deviates somewhat from the observed orientation. 

More recently, Tonry et al. (2000) used the results from the Surface Bright- 
ness Fluctuation Survey to model the local and large-scale flows. Their best-fit 
model includes two attractors, one of which having a best-fit location coinci- 
dent with the Virgo cluster and the other having a fit location slightly beyond 
the Centaurus cluster (which is commonly referred to as the Great Attractor). 
Assuming an attenuated power-law mass distributions for the two attractors, 
their best-fit model has enclosed mass (in excess of background density) of 
7 x 1O 14 M for Virgo and 9 x 1O 15 M for the Great Attractor within spheres 
centered on the attractors with radii reaching the Local Group. For compar- 
ison, the corresponding masses in our simulation are 7 x 10 14 /iQ 1 M Q for the 
Virgo and 8 x 10 15 /ig 1 M for the Centaurus, in good agreement with the 
observational estimate by Tonry et al. 

We follow the trajectories of the above 11 particles up to the time t = to + Qtn- 
The right panel of Figure 4 indicates that the particles have traversed a co- 
moving distance of ~ 5/iq 1 Mpc towards the direction of the Virgo cluster but 
have not fallen into it. The physical distance to Virgo (which is proportional 
to a(t)) would increase exponentially at late times. Neighboring particles show 
a similar behavior, so this conclusion is not likely to be affected by the ran- 
dom sampling of particles from the original data set described in § 2. By 
t = t + 6t H , the peculiar velocity field at moderate overdensities is heavily 
damped (see § 3). Most of the displacement towards Virgo takes place by 
t = t + 2t H - Note that the Virgo cluster itself has also moved by a few co- 
movme; h 1 Mpc in response to the pull of the Great Attractor. We conclude 
that the Local Group will not be gravitationally bound to the Virgo cluster, 
given the initial conditions of the simulation that are matched to the observed 
galaxy distribution. 



5.2 The Milky Way and the Andromeda Galaxies 

The Local Group is a centrally concentrated system in which half of all 
members lie within ~ 450 kpc of its center-of-mass. Most of the mass is 
concentrated in the Andromeda and the Milky Way subgroups. Unfortu- 
nately, the mass estimate of each subgroup suffers from large uncertainties. 
For example, the mass of the Andromeda subgroup ranges from a value of 
M A = (13.3 ± 1.8) x 10 n M o (Courteau & van den Bergh, 1999) to a lower 
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Fig. 4. Particle distribution in a slab of thickness — 15 < SGZ < 15/iq 1 Mpc pro- 
jected onto the supergalactic XY plane. The left panel shows the projection today, 
and the right panel shows it at t = to + . The Virgo cluster appears near the 
top of the image. Open circles mark particles that are within a radius of 2/iq 1 Mpc 
from the supergalactic center at the present time. They advance towards Virgo by 
a comoving distance of ~ 5/j.q 1 Mpc, but do not fall into it at late times. Most of 
this displacement is completed by t = to + 2tff. 

estimate of 7.0±3° 5 5 x 10 n M Q (Evans et al, 2000). The mass of the Milky Way 
ranges from M MW = (8.6 ±4.0) x 10 n M o (Zaritsky, 1999) to a larger value of 
19l?y x 10 n M o (Wilkinson & Evans, 1999). All of the existing estimates are 
consistent with each other to within the quoted errors. Under most estimates, 
the total mass of the Local Group is close to Mlg = (23 ± 6) x 1O 11 M (van 
den Bergh, 1999). 

We may evaluate the present-day overdensity of the Local Group, 5lg, using its 
total mass and the distance of I = 740±40kpc (Binney & Merrifield, 1998) be- 
tween the Milky Way and Andromeda, yielding l + 5 LG = 3M LG /[47r(£/2) 3 p] = 
130/iq 2 ■ The present overdensity of the Local Group clearly exceeds the critical 
collapse threshold of S c = 17.6. 

Adopting the measured radial component of the relative velocity between 
the Andromeda and the Milky Way galaxies v r = —120 km s _1 (Binney & 
Tremaine, 1987) and ignoring its tangential component, it can easily be shown 
that the total energy of the Local Group, ^Hlg v2 ~ GM\M MW /£, is negative, 
where /zlg — Mj^M MW / (M^ + M M w) is the reduced mass of the Andromeda 
and the Milky Way subgroups. This result is obtained for either of the latter or 
the former pair of mass values quoted above for Ma and Mmw- The transverse 
component of the relative velocity is commonly assumed to be much smaller 
than the radial velocity (Einasto & Lynden-Bell, 1982). Based on the radial 
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separation and the relative velocity between the two galaxies, one finds that 
they are likely to merge within a Hubble time. A pedagogical discussion on 
the dynamical evolution of the Local Group based on the least action princi- 
ple can be found in Peebles (1993). Dynamical evaporation of stars out of the 
merger product can be neglected on the time-scale of interest here (Binney & 
Tremaine, 1987). 



6 Conclusions & Discussions 

We have simulated the future evolution of our cosmic neighborhood in a uni- 
verse dominated by a cosmological constant using an N-body code and the 
initial conditions that were reconstructed from the observed galaxy distribu- 
tion of the IRAS 1.2 Jy survey. We find that the large-scale structure and the 
mass distribution of bound objects will freeze in ~ 30 billion years from today 
due to the accelerated expansion of the Universe. The Local Group of galaxies 
will get somewhat closer to the Virgo cluster of galaxies in comoving coor- 
dinates, but will rapidly recede from Virgo in physical coordinates. However, 
the overdensity inside the Local Group is well above the required threshold 
(5 C = 17.6) for it to resist the repulsive gravitational force of the cosmological 
constant. Therefore the Milky Way and the Andromeda galaxies are likely to 
merge within a Hubble time. 

If the Universe is dominated by a cosmological constant, then the Local Group 
is detached from the rest of the Universe and the physical distance from us 
to all other systems that are not bound to the Local Group will increase 
exponentially with time in the distant future. Combining our simulations with 
the results of Loeb (2002), we predict that when the age of the Universe will 
be ~ 100 billion years, there will only be one massive galaxy for us to observe, 
namely the merger product of the Andromeda and the Milky Way galaxies. All 
other systems that are not gravitationally bound to the Local Group will exit 
through our event horizon and their images will fade rapidly while remaining 
frozen on the time of their exit. 

The precision of our quantitative results is limited by the accuracy of the initial 
conditions of the simulation at z = compared to the true mass density field 
in the nearby universe. The locations of the simulated clusters of galaxies are 
shifted by a few megaparsec relative to the recent observational estimates, and 
in principle the match between the two can be improved in the future by refin- 
ing the reconstruction method of the initial conditions of the simulation from 
the observational data. While we believe that the general predictions derived 
in this paper are independent of the limited mass and spatial resolution of our 
simulations, it would be appropriate to perform in the future higher-resolution 
simulations with improved initial conditions that match the observed galaxy 
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distribution. Since the observational uncertainties in the distance measure- 
ments are still substantial, the precision of the observed galaxy distribution 
will also get better with future observations. The combination of both of the 
above improvements will allow us to make more accurate predictions for the 
future evolution of the nearby large-scale structure. 
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